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Abstract 

Employing the currently discussed notion of pseudo-hermiticity, we 
define a pseudo-unitary group. Further, we develop a random matrix 
theory which is invariant under such a group and call this ensemble 
of pseudo-Hermitian random matrices as the pseudo-unitary ensem- 
ble. We obtain exact results for the nearest-neighbour level spacing 
distribution for (2 x 2) -matrices which has a novel form, s log ^ near 
zero spacing. This shows a level repulsion in marked distinction with 
an algebraic form in the Wigner surmise. We believe that this 
paves way for a description of varied phenomena in two-dimensional 
statistical mechanics, quantum chromodynamics, and so on. 



PACS Nos : 05.45. +b, 03.65.Ge 

Postulates of quantum theory require the observables to be represented 
by Hermitian operators as only real eigenvalues correspond to measurements. 
However, it has recently been emphasized that there are certain Hamilto- 
nians describing the quantum systems which possess real eigenvalues even 
though they are not Hermitian. Many of these systems are invariant under 
space-time reflection, i.e. invaiant under a joint action of parity {V) and 
time-reversal (T) [|^, 0, In this context, the concept of pseudo-hermiticity 
was introduced where it was shown that "PT-symmetry is a special case 
of pseudo-hermiticity. Pseudo-hermiticity of an operator or a matrix O is 
simply defined through the condition : 0^=r}Or}~^ with rj a metric and ^ 
representing the usual adjoint or conjugate-transpose. Remarkably, it was 
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subsequently shown that non-VT invariant systems that possess real eigen- 
values are also pseudo-Hermitian Physical situations of great interest 
belong to the above discussion. This includes two-dimensional statistical me- 
chanics where parity and time- reversal are broken (preserving VT) p, |[, 
quantum chromodynamics where chiral ensembles are used to describe the 
statistical properties of lattice Dirac operator , spin-rotation coupling lead- 
ing to an anomalous g- value for muon , and related fields. In this Letter, 



we present a random matrix theory which describes spectral fluctuations in 
systems which are pseudo-Hermitian and pseudo-unitarily invariant. The 
two aspects which are particularly notable are the simplicity of this novel 
description and the fact that this theory is natural when parity or (and) 
time-reversal is (are) violated. 

The problem of two-dimensional statistical mechanics is obviously con- 
nected with anyon physics and hence to the behavior of an electron in an 
Aharonov-Bohm medium i.e. a medium filled with non-quantized mag- 
netic fluxes, reminiscent of the theory of fractional quantum Hall effect |12 



Important to note here is also another motivation which stems from a specu- 
lation due to Nambu that this might serve as a model for theoretical ideas like 



the quark confinement in a medium of monopoles |T3[. In this context, it is 



known that the spectral fiuctuations of an Aharonov-Bohm billiard exhibits 
an interpolating behavior with respect to the strength of the fiux line [in . 



These billiards are experimentally realized in terms of quantum dots in the 
presence of fiux lines. It is of great interest to find an appropriate random ma- 
trix description for such PT-invariant systems. Pseudo-hermiticity appears 
in several contexts. It is instructive to note that in the mean-field, RPA 
description of nuclei [1^, the stability matrix leading to an eigenvalue prob- 



lem can be checked to be pseudo-unitary. In the context of regularization of 
quantum field theories, pseudo-hermiticity and the associated improper met- 
ric was used by Dirac [T^, Pauli |T^, and particularly by Gupta and Bleuler 
jTB], and others |T^. Let us first establish the pseudo- unitary symmetry. 

Consider vectors x and y residing in a vector space V and a fixed metric 
rj. In this vector space, we define a pseudo- inner product (57- norm), which 
can be written in the usual quantum mechanical notation as {x.\r)'y). We shall 
consider symmetry transformations which preserve the J^-norm between the 
vectors. We consider the Cayley form, D = e*^ as a symmetry transformation 
acting on x, y where G is pseudo-Hermitian in accordance with r]Gr]~^ = G^ 
Noting an interesting feature of D : 

= rje-'^rj^ = riD^ri~\ (1) 
let us call D as pseudo-unitary with respect to rj. rj equal to unity makes 
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D unitary trivially. To establish that D is indeed a symmetry transforma- 
tion, we need to show that the transformation preserves the 77-norm and a 
consistently-defined matrix element. 

Let us assume that x (y) — >• x' (y') = Dx (Dy). Then, the pseudo- unitary 
symmetry is defined by preserving the pseudo-norm : 



(x'lw') = (Dx|7/Dy) = (x|7/y). (2) 

In proving (2), use e'^^'^'qe^^ = e~**^^r/e*'^r/~^r/ = e~**^^ 6**^^77 = rj. Under 
the same pseudo-unitary transformation, the matrix element of an arbitrary 
operator. A, transforms as 



(x'lT^A'ly') = (x|77A|y) 



(3) 



if D A D"^ = A'. 

Let us now prove that pseudo-unitary matrices form a group under matrix 
multiplication. For closure, let Di and D2 be two pseudo-unitary matrices. 
D1D2 is pseudo-unitary because r] 
It easily follows that D 



-1/ 



DiDa)^?/ = r] f/T/ ^D\r] = (D1D2 



-1 



unitary : r) (e" 



V 



^ is pseudo- unitary with respect to 7/ if D is pseudo- 
vj jj _ r^Yie identity matrix acts as the unit 
element of the symmetry transformation. Finally, since the associativity is 
guaranteed, the N x N pseudo-unitary matrices form a pseudo-unitary group 
of order N, PU{N). 

In the following, to keep the proceedings simple and explicit, we consider 
Hamiltonians in their matrix representations. Also, in the spirit of the origi- 
nal work of Wigner we consider (2 x 2) matrices as they bring out most 
of the essence. In this context, there i a recent generalization of Wigner sur- 
mise for 2x2 matrices |^T|. Thus, we concentrate on PU{2) and consider 
the following pseudo-Hermitian matrix. 



H = {H,,} 



a 
ic 



-ib 



(4) 



a, b, c being real. Consequently, e**^ will be a pseudo-unitary matrix. For the 
above matrices, a metric is 



-1 

1 



(5) 



This metric may be interpreted as the parity operator V, and the complex 
conjugation, JCq as time-reversal operator T. With these operations, it may 
be verified that H is PT-invariant in addition to being P-pseudo-Hermitian. 
Besides these commuting V and T operators, if we choose T as the Pauli 
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matrix times the complex conjugation, Ux/Co, they do not commute, however 
preserving other conclusions. 

This group admits three generators and an identity, viz.. 



Pi 
Pi 



1 
i -1 

-1 
1 



P2 
,1 



1 -i 

-1 

1 
1 



(6) 



Note that H = al + hpi + cp2 - {h + c)p3. It is interesting to see that pi and 
P2 are pseudo-hermitian and pseudo- unitary, possessing eigenvalues ±1. It 
may be recalled that the Pauli matrices cr^ and ay are Hermitian and unitary. 
Further, the generators satisfy the following important properties. 



pI 
\pi,Pj 



pI 



pI 



I, 



(7) 



with = = CI3 = CI3 = Cg^ = C|i = 2 and Cf^ = 5. All the structure 
constants can be found with the help of commutation relations and symmetry 
properties, and they turn out to be ±5, ±2, or 0. Interestingly, the following 
relations between the structure constants hold : 



^kl — ~^lk 



(8) 



m=l 



thus making it a Lie group and defining a Lie algebra p3 



We now consider a Hamiltonian H which is diagonalizable by D, i.e.. 



H 



E+ 
E_ 



.D 



(9) 



The eigenvalues of H are a ± ybc {be > 0). The corresponding matrix, D, 



1 i/r 
ir 1 

is pseudo-unitary under the metric, 

" 1 
^=10 



(10) 



'11^ 
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The eigenvalues are 



E4 



a± 



c br 
2^ ^ 2" 



(12) 



where r = \/c/b (0 < r < 00). 

Consider that the matrix H is drawn from an ensemble of random matrices 
with a Gaussian distribution given by |2( 



P(H) 



_ 1 tr H^H 



Accordingly, the joint probability distribution of a,b,c is 



P(a, 6, c) 



1 



-e 2<T 



2(7ra2)i 

From (4) and (9), we have the following relations 



(13) 



(14) 



2^ 



r{E+ - E^ 



(15) 



The Jacobian, J connecting {a,b,c) and {E^,E_,r) is ^^'^2r^ ^ ■ With these, 
the joint probability distribution function (j.p.d.f.) of eigenvalues is 



P{E+,E^ 



\E+ - E_\ [ {E+ - E_] 



2f7ra2)^ 



4a2 



7:^ — 



(16) 



Following the Dyson Coulomb gas analogy, this j.p.d.f. can be written as an 
equilibrium distribution of two interacting particles with a partition function 
P{E^, E_) 2(xi, X2) = e~^^*^^'i'^2)_ ii is interesting to note that the TC has 
a potential term involving the logarithm of the modified Bessel function along 
with the familiar harmonic confinement and the two-dimensional Coulomb 
potential. 4cr^ plays the role of inverse scaled temperature. 

Integrating with respect to E_ gives the average density, shown in Fig. 
1. This is not amenable to an analytically closed form. 

Perhaps the most well-studied characterizer is the nearest-neighbour level 
spacing distribution, P{S). This gives the frequency with which a certain 
spacing between adjacent levels occurs. For the Wigner-Dyson ensembles, 
P{S) ~ S^°e~^^ where (3q is 1, 2, and 4 for the orthogonal, unitary, and 
symplectic ensembles. A wide variety of systems display universal proper- 
ties possessed by random matrix ensembles as can be seen in |20, |23, |25 



However, there are systems that display intermediate statistics [^, ^ ^ . 
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These systems range from examples of billiards in polygonal enclosures, three- 
dimensional Anderson model at the metal- insulator transition point, and so 
on. On the other hand, there have been important developments on non- 
Hermitian ensembles since long where the eigenvalues are complex pO|, p5| , 
and where an ensemble of unstable states are considered [^. Clearly, the 
ensemble developed here does not fall into any of the known categories and, 
indeed, displays some novel features as shown below. 

The spacing distribution, P{S), is given in terms of the j.p.d.f. by 

/oo roo 
/ P{E+,E_)5{S -\E+- E_\)dE+dE_ 
-oo J —oo 



TTCr 



This result is distincty different and very interesting (Fig. 2), particularly 
for its behaviour near zero spacing. Near 5 = 0, the probability distribution 
varies as S'log ^. This follows from the asymptotic properties of the modified 
Bessel function. 

We present the form of two-time correlation function for a complex sys- 
tem with spectral properties given by Gaussian pseudo-unitary ensembles 
(GPUE). In this context, we consider a system with a Hamiltonian H G 
GPUE, and an observable given by an operator, V G an another GPUE. 
Imagine the system to be in thermodynamic equilibrium with a canonical 
density matrix, p = e~^^. Following an extensive study along the lines of 
pO| , the time correlation function, C(t) = Z~^{P) tr pV(0)v(t) decays over 
long times ~ (^!£ii/£)3_ 

Finally, we wish to point out an aspect of general importance, encoun- 
tered on many occasions in many-body theory. To give one concrete example, 
in the theory of collective excitations of Fermionic systems, a mean-field de- 
scription is used where a collective state is first expressed in terms of particle- 



hole excitations [^. Here, one generally encounters a matrix equation like 
H\E' = A$ with H a Hermitian or unitary operator. The above problem may 
be transformed into an eigenvalue problem for H', i.e., H'$ = A$ with H' 
a pseudo-Hermitian or pseudo-unitary operator. With this, there are many 
results immediately possible. First of all, the eigenvalues will either be real, 
complex-conjugate-pairs, unimodular, or they occur in pairs such that prod- 



uct of eigenvalues is unimodular ||2^. Secondly, the statistical properties 



of the eigenvalues related to collective excitations will be distributed in ac- 
cordance with the results obtained for GPUE above. Thirdly, there will be 
long-time tail for relaxation in the same spirit as the time correlation. 

The above results are found for 2x2 matrices. Although for N x N 
matrices, invariant under PU{N), the results are not known, we conjecture 
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that the fluctuation properties will have a similar form as above. As dis- 
cussed earlier, the result found here gives a new universality corresponding 
to systems which are pseudo-unitarily invariant. In such systems, parity and 
time-reversal may be individually broken, preserving their joint action. This 
universality also includes those pseudo-Hermitian quantum systems where 
VT is broken. The examples discussed include quantum chromodynamics, 
two-dimensional statistical mechanics, and so on. 
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FIGURE CAPTIONS 



1. The average level density of an ensemble of 2 x 2 random Gaussian 
pseudo-unitary ensemble is shown here. 

2. The nearest-neighbor level-spacing distribution is shown here. For 
comparison, the results corresponding to the Wigncr-Dyson ensembles 
corresponding to orthogonal and unitary symmetries are also shown. 
Whereas the level repulsion is hnear and quadratic in the orthogonal 
and unitary ensembles, here it is of the form slog(l/s), as shown in 
the inset. This then is a new universality. 
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Fig. 2 



